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Abstract 

We study the propagation of electromagnetic (EM) waves in a slab-type waveguide which walls 
consist of three-dimensional topological insulator (3D TI). The results show that a cutoff frequency 
with topological stability limits the spectrum that propagates along the waveguide and are in 
agreement with experimental observations. Our approach also provides a way to measure the 
penetration length of surface metallic states in 3D TI. 
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I. INTRODUCTION 


The spin quantum Hall (SQH) electronic state was discovered experimentally in the 
HgTe/CdTe quantum wells being called two-dimensional topological insulator (2D TI) PJ- 
Sj. Subsequently, this electronic phase was found in materials such as Bi 2 Tes and Bi 2 Se 3 , 
named 3D TI [J, 5j. They are characterized by a full insulating gap in the bulk and gapless 
metallic surface (edge for 2D TI) states which are protected by time reversal symmetry (TRS) 
J6j. These surface states remain robust against small non-magnetic perturbations and they 
are characterized by topological invariants that are preserved as long as hamiltonian varies 
smoothly [7]. In addition to their fundamental interest, their special properties make these 
states useful for several applications, from spintronics to quantum computation and optics 

mm- 

The long wavelength regime of EM wave of 3D TI may be described by a topological 
field theory (TFT) whose coefficients are quantized in terms of the fine structure constant, 
a = e 2 /hc |12j . This TFT describes the topological magnetoelectric effect (TME) on the 
surface of the material when the TRS is broken at the surface [T3] . Thus, in the presence 
of EM fields, an electric field induces a magnetization whereas a magnetic field induce an 
electric polarization. This topological effect is caused by a quantized Hall current induced 
on the surface of the TI in response to an electric field, in (3+1) D, the TFT is given by the 
action below [1%]: 

s= iS {tp ~m ,i ' x+ ^s e ' sdix - w 

where E and B are the electric and magnetic fields respectively, /i and e are the permeability 
and permissivity of the TI and 9 is the relevant topological parameter accounting for the 
TME effect and assumes 9 = 0 for an ordinary insulator and 9 = n for TI. This effective 
theory is valid only if a surface gap is open on the TI, this can be done by covering its surface 
with a thin magnetic film or applying an external magnetic field perpendicular |T6j. Among 
other consequences, when light is send to a TI border its behaviour is deeply modified, for 
example, the light polarization plane is rotated by a universal angle leading to changes in 
the usual Kerr and Faraday effects [13] . 
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II. SLAB WAVEGUIDE WITH TOPOLOGICAL INSULATOR WALLS 


Here, we study the propagation of EM radiation confined to a slab-type waveguide which 
walls are the surface borders of a 3D TI. The interior of the guide is filled with an usual 
dielectric (vacuum, for simplicity), with this, there is a variation in 6 at the interface of 
the TI and thus a contribution from topological term in Eq. ([Tj) . We consider an aperture 
L along x and an infinite extension along y and z directions. In addiction, an external 
magnetic held is applied perpendicularly to the walls, as schematically drawn in Fig|l]. 



FIG. 1: Slab waveguide with aperture length L and TI walls. In 0 < x < L the slab is filled with an 
ordinary dielectric, k! is the wave vector of a mode allowed by the system, Ju is the Hall current 
induced by EM waves and circulates along the interface vacuum/TI. h is an unitary vector normal 
to TI surface and gives the direction of the external magnetic field applied (Bo)- 

The functional variation of Eq. ([!]) yields Maxwell’s equations in the presence of TME, 
that read: 


V ■ E = — 
e 

( 2 ) 

J j 

V- B = 0, 

( 3 ) 

1 - 

V x E = — d t B, 
c 

( 4 ) 

V xB = fiJ H + —d t E. 

( 5 ) 
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Most studies on TME consider the charge (pn) and Hall current (Jh) densities perfectly con¬ 
centrated on TI surface using a Dirac 8(f) distribution. We take into account the penetration 
length, l, of surface metallic states on TI bulk. Thus, in our approach: pn = — n-B ) and 
J H = yf (n x E), where l ~ 10 _9 m (few nanometers) for usual 3D TI, like Bi 2 Se 3 [TT| H8 j. 
While fi is a unitary vector perpendicular to the walls of the waveguide. 


III. RESULTS AND DISCUSSION 

Once this slab waveguide conhnes waves only in x direction, we assume the modes are 
propagating freely along z, so it is necessary to solve the Maxwell equations for planar 
solutions of EM fields such: E(f ., t) = E(x)eJ kz ~ u 9 and B(f , t) = B (x)eJ kz ~ ut \ Taking then 
into Eqs. Q and ([5]) we obtain the transverse fields: 

E x = 2 l _ {c 2 kd x E z + cwdyB z + pcuj9'(n x E)^j , (6) 

By = - y~ —W 2 ( c2 kd y E z - cud x B z + pcu9'(h x E) y ) , (7) 

fJjCUJ C rC \ / 

B x = l _ ^2 (c 2 kd x B z - peaud y E z + pc 2 k9'(h x E)^j , (8) 

By = 2 l _ ^ 2 (c 2 kd x E z + peaxdyBz + pc 2 k9'(n x E)^j , (9) 

where 9' = ^ and n — x (see Fig. [TJ) . For the slab waveguide is only necessary to analyze 
the dynamics over z, so, from Eqs. ([2]) and ^ we obtain: 

d 2 x E z (x) + uj‘ 2 - k 2 ) E z (x) - & -d x B z (x) - p9 ,2 E z (x) = 0, (10) 


d 2 B z {x ) + - fc2 ) B z (x) + p9'd x E z (x) = 0. 


( 11 ) 


To decouple these equations we must depart to higher orders. Doing that for B z (x), we get: 


pe 


9 , 


pe 


dt + 2 ( ^u 2 -k 2 - '-9 rz )d 2 +(^u 2 -k 2 )[ ^u 2 ^ f 0 

1 e kc 


pe 


p id 


1/2 


B z (x) = 0, (12) 


where E z (x) is obtained from Eq. Taking 9 = 0, we recover the standard equations as 
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expected. The solution has the form: B z (x) = Bie~ 1+X + B 2 e~ 1 ~ x , valid for the x > L, for 
x < 0 the exponential signals are inverted so that solution remains finite. The constants B\ 
and B 2 are fixed field amplitudes within the TI, while: 


■ 2 pe 2 l 1 qi 2 _j_ \Jpc 5 9 ,2 k(2ec 3 k 3 + pc 3 9’ 2 k — ec 2 k 2 ui — 2 pe 2 cku 2 + pe 2 oj 3 ) 

^ c 2 ^ e ec 4 k 


1/2 


(13) 


From the wave equation for EM fields, j± must satisfy 7 ^+ 7 \ — k 2 -P ^oj 2 = 0, so that we 
finally obtain: 


kc 

~ = 
to 


1 - ( — 

OJ 


2 2 c cud 

, u e = - r . 

7r e l 


(14) 


Waves with cn < do not propagate in TI surface, since k becomes imaginary resulting in 
attenuation of the EM radiation. These frequencies propagate along z direction inside the 
waveguide only, therefore, uig is the topological cutoff frequency, such behaviour is illustrated 
in Fig. [ 2 ] for Bi 2 Se 3 . 



FIG. 2: Normalized wave vector in arbitrary units (a.u.) as function of u. We consider l = 2 
nm and Bi 2 Se 3 dielectric parameters, p = 1 and e = 100 m, so the cutoff frequency is ug « 10 13 
Hz. 


Since the cutoff frequency depends on /, the system proposed here can be useful to es¬ 
timate the penetration length of the surface states in 3D TI bulk by measuring ujg. The 
topological nature of utg stems from the TME on the surface of TI and this frequency corre¬ 
sponds to the minimum energy required for EM radiation through the surface of TI. Thus, 
EM waves with uj < ojg are reflected keeping pn and Jh on the slab surface and propagate 
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along the waveguide. These results are in agreement with experimental data reported in 
Ref.[2D], which were made measurements of the reflectivity and optical conductivity of four 
3D TI (including Bi 2 Se 3 ) from 5 to 300/1 and from sub-THz to visible frequencies. The 
data show a sharp drop in reflectivity (R(uj)) at 00 ~ 10 13 Hz and R(u>) ~ 1 for u> < c jg, as 
provided in our approach. 

Along x, the wave vector in Tl can be expressed as a function of cn and sweeping the 


(15) 


same frequency values made for Eq. (14), we found that Eq. (13) takes the form: 7 ±(cu) = 
£((n) ± ir{uj). Therefore, the wave vector perpendicular to Tl surface have real (£) and 
imaginary (r) values. Typical behaviours of £(cu) and t(u) are shown in Figs. [3] (a) and [3] 
(b) respectively. Below the cutoff frequency, the EM wave is strongly attenuated when it 
penetrates the TI, once £(cj) 7 ^ 0, so when u: < loq there is an efficient confinement of the 
radiation along x, according to nontrivial boundary conditions on EM fields (obtained from 
[2] and [5]) : 

n • ( E vac — eE T i) = p H 
n x (B vac - 2 xl) = J H , 

where subscript vac denote the EM fields in vacuum. I 11 addiction, the stationary wave 
condition require that inside the wavequide, the wave vector along x (k x ) is discretized and 
depends on the geometrical parameters of the waveguide, i.e., k x = nn/L (n = 1,2,3,...). 
Due to the large attenuation of EM radiation around cog, the skin depth (l/£) of EM wave 
in TI decreases rapidly in u 0 ~ ujq. For u> > t jg, vanishes, while t(oj) reverse signal 

abruptly resulting in transmission of EM radiation through the walls. This is in accordance 
whit optical conductivity (cr(o;)) measurements performed in Ref.p0], where it was observed 
large values for cr(u) in this frequency regime. 

For the sake of completeness, let us give a brief analysis of transverse polarization modes. 
The transverse magnetic (TM) mode, B z (x) = 0, so, from Eqs. (10) and ( [TT| we get: 
d 2 f E z (x) + ((fcj 2 — k 2 — p#' 2 ) E z (x) = 0 and d x E z (x ) = 0. Adopting procedure previously 
applied, we have cvg ~ 10 13 Hz as EM waves with generic polarization (B z (x) 7 ^ 0 and 
E z (x) 7 ^ 0). In the transverse electric (TE) polarization, E z (x) = 0, then: d x B z (x ) + 
((^cu 2 — k 2 ) B z {x) = 0 and d x B z (x ) = 0. Note that for waveguide operating in the TE mode 
there is no contribution of topological parameter 9, the transverse components (Eqs. 


of E and B fields vanish on the walls and inside the TI leading to pn = -^j(x-B) = —^-B x = 
0, that is, TE waves do not produce TME in the slab waveguide. With theses conditions, 
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(b) 
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FIG. 3: Evanescent and oscillatory components of the perpendicular wave vector, 7 +(w) = £(ui) + 
We take l = 2 nm then u ig « 10 13 Hz. When uj = uig, (a) £(w) vanishes abruptly and (b) 
for t(u>) there is an abrupt reversal of signal. 

the cutoff frequencies along the slab are similar to a waveguide with perfectly conducting 
walls, i.e., oj m = mirc/L (m = 1,2,3,...). For EM waves in the TE mode, frequencies less 
than uj m do not propagate along the slab waveguide similarly to the conventional waveguides 
described in Refs. mm- 


IV. CONCLUSIONS 

In summary, we find that waveguide made from 3D TI walls can propagate EM radiation 
efficiently and have cutoff frequencies with topological stability, which value depends on l of 
the metallic surface states. The value provided here for uj e agrees with experimental results 
of optical reflectivity and conductivity of 3D TI reported in Ref. [120]. This slab waveguide can 
be useful as a prototype to determine a microscopic quantity, l, by measuring a macroscopic 
one, ojg. In the TM mode, the cutoff frequency has the same characteristics and value 
of the generic polarization. On the other hand, TE mode has no topological properties 
due to TME not manifest for this polarization state, so the cutoff frequency depends only 
on the geometric parameter (L) similar to waveguide with usual conductive walls. Other 
geometries as cylindrical waveguide, coaxial cable, resonant cavities, etc. can reveal new 
and useful electrodynamics properties of 3D TI. 
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